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Context. In the spectra of variable pulsating stars, especially Mira stars, the detection of intense hydrogen emission lines has been 
explained by the presence of a radiative and hypersonic shock wave, periodically propagating throughout the stellar atmosphere. 
Previous observation of the Mira star omicron Ceti around one of its brightest maximum light led to the detection of a strong level of 
^\ ' linear polarization associated to Balmer emissions, although the origin of this phenomenon is not fully explained yet. 

^SJ , Aims. With the help of spectropolarimetry, we propose to investigate the nature of shock waves propagating throughout the stellar 

atmosphere and present, for omicron Ceti (the prototype of Mira stars), a full observational study of hydrogen emission lines formed 
in the radiative region of such a shock. 
' Methods. Using the instrument NARVAL mounted on the Telescope Bernard Lyot (TBL) in Pic du Midi Observatory (France), 

, we performed a spectropolarimetric monitoring of omicron Ceti during three consecutive pulsation cycles. For this survey, the four 

Stokes parameters (I for intensity, Q and U for linear polarization, and V for circular polarization) were systematically collected, with 
a particular emphasis on the maxima of luminosity, i.e. when a radiative shock wave is supposed to emerge from the photosphere and 
starts to propagate outward. 

Q , Results. On hydrogen Balmer lines, over a large part of the luminosity cycle, we report clear detection of polarimetric structures in Q 

5_j ■ and U Stokes spectra (and also in V Stokes spectra but to a lesser extent). We report a temporal evolution of these spectropolarimetric 

signatures, which appear strongly correlated to the presence of an intense shock wave responsible for the hydrogen emission lines. 
' We establish that the hydrogen lines are polarized by a physical process inherent to the mechanism responsible for the emission line 

formation: the shock wave itself. Two mechanisms are thus considered: a global one that implies a polarization induced by some giant 
convective cells located around the photosphere and a local one that implies a charge separation due to the passage of the shock wave, 
inducing an electrical current. Combined with the existing turbulence, this may generate a magnetic field, hence polarization. 
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Received July 2011 and t hermodynamic couplin gs between convection and oscilla- 

. tions. IXiong & Dend ( l2007h found that turbulent pressure plays 

' a key role for the excitation mechanism. The typical period of 

1 . Introduction Mira stars' pulsation extends from 150 to 1 000 days, the period 

during which the luminosity might change up to 2.5 visual mag- 

Mira stars are cool, evolved and variable stars. Its prototype nitudes on average. It has to be noted that those stars are not 

is omicron Ceti (as known as Mira), which is the object of this as regular as, e.g., RR Lyrte. Therefore the duration of a cycle 
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paper. On a Hertzsprung-Russell diagram, Mira stars are associ- might vary with time 

lJ ated to the asymptotic giant branch (AGB). Indeed, their main- r^, , , , ^ t^n , 

^ • 1 . 1 u- u 11 »u The low surface temperature of those stars (~3 000K) al- 

sequence mass is less than ten solar masses, which allows them , ^ , . ^ , , • ■ ■ , 

' ' ' u \r^n u u u j j u i- u • lows for the cxistcncc of molecules in the very extended at- 

to reach the AGB phase where hydrogen and helium burn in , , • , • , 

, „ , f u J T-u J mosphere and in the circumstellar envelope, all of which are 

shells around a core of carbon and oxygen. They undergo ra- . , r r , , ,- . 

J. , , . ... 1 • J u 1 • typical features of late-type stars. For the Mira stars of the M 

dial pulsations whose origin is explained by several theories, , , , , . , 

^, . , . . , . lhno< SK u- u spectral type, the most present molecules are oxygen-rich ones 

the main one being a supposed ;<--mechanism OCoxll 19801) . which ^ , r^.^ ™ , , , . . r 

, , , 1 f • 1 1- • c Tc J 1 such as H9O, TiO, or CO. These molecules are at the origin of 

leads to a cycle of visual light variations. Some self-excited pul- , " , , , • , ,,,,,,, 

, , fT,,. 1 J t/T"! — A t 1 I hnnol A e the Strong molecular absorption undergone by th e black-body 
sation models of Miras already exist (llreland et al.ll2008l and ret- ° ,. ^ .. . i r,.' ■ i Nr.-ri I I n^ ■ 1 

*u • N u » »u u • ui f 1 spectrum at this temperature .Tsuii. .19oa, Tsuij 1971al Tsuij 

erences therein), but the mechanism responsible for the pulsation i^'lJZ., 1 rr 1 ,-.^'; . i 7J0r.r.L r.^. . — , — : — . T 

, A . f^ A . f • . . 1971bi J0rgensedl994 and Gilletll988h . This absorption is such 

has not been identmed yet. The free input parameters are mass, '-; f^r T , T r , 

, . -^^ • ^ u 1 ^. 1 • • 1 *u that the observed spectrum departs a lot from the correspond- 

luminosity, composition, microturbulent velocity, mixing length ■ , • ^ , 

-.IT- 1 w J J . ing black-body spectrum. Furthermore, and especially for Mira 

and turbulent viscosity. Using a non-local time-dependent the- ° -' ,^ , , , , ' X , , , 

, , , \. V- T~^\^r:no^ u J . stars. Strong emission has been obscrvcd on Balmer hydrogen 

orv of turbulent convection. Xio ng et all (119 98) showed that the ,. , . ° , , ,. , , . . . , 

, ^. , . u i-. f iv/r- » ■ u J • lines, lasting during up to 80% of the luminosity period, 

pulsational stability of Mira stars is governed by the dynamic > & & r j t- 

In past decades, several studies - based on high resolution 



* Based on spectropolarimetric observations obtained at the spectroscopic surveys - have focussed on the atmospheric dy- 
Telescope Bernard Lyot (TBL at Observatoire du Pic du Midi, CNRS namics of several classes of radially pulsating stars located along 
and Universite de Toulouse, France). the instability strip: from the population II cepheid W Vir stars 
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("Lebre & Gille3ll992h up to the RV Tauri stars ("Lebre & Gillet" 
1991), the Mira stars |Mvai-ez et al. 2000 and Alvarez et al. 
20011) . and the post-AGB stars (Lebr e et al.lll996l) . It is now ad- 



mitted that radially pulsating variable stars host radiative hy- 
personic shock waves that propagate periodically (the period is 
hardly changing with time) throughout the stellar atmosphere. 

Two main properties of these shock waves have been investi- 
gated: (i) their radiative nature, generating e mission lines forma - 
tion process (hydrog en and helium lines, see lGillet et alj|T983l ; 
fluorescent lines, see IWillsonI 19761) . and (ii) their hypersonic na- 
ture, from an estimation of their propagation velocity that always 
appears higher than the sound velocity in such media (several 
Mach numbers ; see iGillet et anil998al and re ferences therein). 
Their structure has also been investigated (seel Fadevev & Gilletl 
|2004| and references therein), and three main areas constitute 
those radiative shock waves. First of all, the radiative precursor 
is the area before the shock where the matter has not yet been af- 
fected by the raising wave compression. Then comes the shock 
front. In a length of only several mean free paths, a part of the 
kinetic energy is brutally converted into thermal energy, which 
creates strong temperature, pressure, and velocity gradients. At 
the end, the radiative wake is located behind the front. There, 
the partially ionized matter gets cooler and recombines. The in- 
tense emission lines observed in the spectra of radially pulsating 
variable stars, such as M ira stars, are supp osed to be formed just 
behind the shock front (IGillet et alJI 19831) in the so-called radia- 
tive wake of the shock. 

These studies have also related the shock wave propagation 
phenomenon in the atmospheres of variable stars to their pulsat- 
ing process (which is the /f-mechanism for the large majority of 
these objects) or to more ch aotic behaviour (fo r RV Tauri and 
Post-AGB stars, for example. iBarthes et alJl200 0). These works 
have mainly established that shock waves, which are connected 
to the pulsation mechanism, are emerging from the photosphere 
once a cycle, around the maximum light, and are thus propagat- 
ing periodically throughout stellar atmosphere, ruling its dynam- 
ics, and extending it to several tens of solar radii, progressively 
generating, and shaping a circumstellar envelope around coolest 
objects. 

However, to date, the magnetic nature of these shock waves 
has never been identified or deeply investigated, while since the 
60s, theoretical studies of shock waves propagating throughout 
gaseous media of atmospheric density predict there is an elec- 
tric field (therefo re a magneti c field) in the post-shock region 
(see for instance: lJaffTin|[l965l : iLu & Huang|[l974l : et alj 
Il995h . Moreover, because of the complexity of the problem, this 
magnetic nature has never been introduced in shock structure or 
shock propagation modelling. Also, its impact on astrophysical 
processes, such as radial pulsation or chaotic behaviour, atmo- 
spheric dynamics, circumstellar envelope formation and struc- 
ture, mass loss, and grain formation, has never been investigated 
so far. 

In the late 70s, narrow- and broad-band polarimetric ob- 
servations were performed on cool evolved stars hosting cir- 
cumstellar envelopes in order to investigate the grain distribu- 
tion in stellar environments. With this technique, in a study de- 
voted to the one-year-p eriod variable Mira star, omicron Ceti, 
iMcLean & Covnel([T978 ) were the first to detect a strongly linear 
polarized flux in the emission lines of the Balmer series. Based 
on just a single observation limited to the maximum luminos- 
ity, the polarization in the hydrogen lines (from Hfi to HIO) was 
noticed as two to three times greater than the polarization mea- 
sured in the nearby continumm. To date, this observational fact 
is still an enigma. It has also never been reproduced on any Mira 



star or on another pulsating star, or even characterized during 
a complete pulsating cycle. Moreover, recent modelling of the 
linear polarization of optical emission in red giants dPadevevI 
2007) which investigates light scattering by dust grains in ax- 
isymmetric circumstellar envelopes, has even failed to explain 
such a high polarization rate. 

Recently, a new generation of instruments has been made 
available: spectropolarimeters. Over a large spectral domain, 
they combine high-resolution spectroscopy and information in 
polarized light, and they open new observational windows for as- 
trophysical objects, mainly focussing on stellar magnetism. With 
these instruments, it is now possible to collect spectral informa- 
tion through the four Stokes parameters (V, characterizing circu- 
lar polarization ; Q and U, characterizing linear polarization ; I, 
the classical spectrum) 

[Harrington & Kuhiil (1200 9) present spectropolarimetric ob- 
servations of Post-AGB stars and RV Tauri stars and report, 
for most of these objects, the detection of linear (Stokes Q and 
Stokes U) polarimetric signatures associated to the Ha lines. 
In most objects, the morphology of the polarization feature 
presents short timescale variability that suggests the presence 
of absorbing polarizing gas close to the star, ac c ording to the 
optical pumping model described by iKuhn et aP (l2007h Stong 
mass loss, progressively shaping circumstellar environments, is 
very common in these stars, and the shock waves propagating 
throughout their at mosphere and in duci ng ballistic rn otions has 
also been reported dLebre & Gilletll 1 99 llllibre et al.ll 19961) . 

It is thus now possible for the very first time to estimate 
the linear and the circular polarization evolution in the emis- 
sion lines of Mira stars. As we are convinced that the hydrogen 
emission lines are polarized by a physical process inherent to the 
mechanism responsible for the emission line formation (i.e, the 
shock wave itself), a spectropolarimetric survey performed dur- 
ing the cycle of pulsation of Mira stars appears to be a good tool 
for tackling this problem. In section 2, we present our observa- 
tions of omicron Ceti as collected with the NARVAL spectropo- 
larimeter. In section 3, we interpret the obtained results in term 
of shock wave characterization. In section 4, we draw a conclu- 
sion from this spectropolarimetric analysis. 



2. Spectropolarimetric monitoring of omicron Ceti 

2.1. Observations 

We performed a spectropolarimetric monitoring along three 
consecutive pulsation cycles of the bright star omicron Ceti, the 
prototype of oxygen-rich Mira stars. During one of its typical cy- 
cles of luminosity (with a mean period of pulsation of 33 1 days), 
the visual magnitude of omicron Ceti presents an average vari- 
ation amplitude of eight and its spectral type varies from M5e 
to M9e (see General Catalogue of Variable Stars and AAVSO 
website). Omicron Ceti has the peculiarity of hosting a compan- 
ion (VZ Ceti) at ~70 A.U. believed to be either a white dwarf 
(|Karovskaetal.lll997l lK arovs kal l2006l) or a mam sequence star 



( Ireland et al T l2007h . We consider it to be too far to have any 



influence on the dynamics of the lower atmosphere of omicron 
Ceti. 

The spectropolarimetric data come from the NARVAL in- 
strument mounted on the Telescope Bernard Lyot (TBL) in the 
Pic du Midi Observatory, France. The polarimeter is made with 
optical elements (a series of half-wave, quarter-wave, and half- 
wave rhombs) that analyse the c ircular or linear polarization 
state of a signal dPonati et al.ll997i) . The processed signal is thus 
transmitted in an achromatic way via an optical fiber to the cross- 
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dispersion spectrometer doing the spectral analysis over a very 
wide part of the spectrum (from 375 to 1 050 nm) and over 40 
orders. The polarimetric mode gives us a spectral resolution of 
68 000, and the velocity bin is then Av ~ 4.4 km. s" ' . 

The monitoring was from the 5 September 2007 to the 10 
February 2010, covering three luminosity cycles of omicron 
Ceti. To provide phases during these cycles, we used as heliocen- 
tric Julian date of reference at luminosity maximum (see AAVSO 
website): HJDi- = 2454 150. The phase (/? is then the ratio be- 
tween HJDo - HJDr and the period (see Fig.[T]and Table[T]). 

Fig. 1: Light curve of omicron Ceti (from AAVSO website) for 
the three cycles of the survey. Each blue cross is one spectropo- 
larimetric observation of our survey. 



Table 1 : Journal of our spectropolarimetric observations of omi- 
cron Ceti with Narval. 



AAVSO DAIA FOH OMl CtT - WWW.AAVSO.ORG 




2454492 2454766 2455040 

Ju[ian Date 

Visual Validated • Visual Prei/alidated 



2.2. Data treatment 

The Libre-ESpRIT software dPonati et al.l 1 1997 1) process- 
ing data collected by NARVAL allows us to obtain polarimet- 
ric spectra. The outputs of Libre-ESpRIT process are the four 
Stokes parameters divided by a pseudo-continuum: 

/ 



- the intensity y , 

- the linear polarization 



Q hr - ho' 



U ^ hy - /i35° 

Ic Ic 



- the circular polarization 

^ V ^ Iq-Iq 

Ic Ic 

(I o : left-handed circular polarization, Iq : right-handed circu- 
lar polarization). 

From this information, we can deduce the rate and the angle 
of linear polarization (see Bagnulo et al. 2009): 
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The Libre-ESpRIT software automatically corrects a 
pseudo-continuum. It deduces a black-body continuum from the 
observed spectrum by which this spectrum is divided. This fea- 
ture is quite a problem in our case because the spectra of cool 
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Mira stars are very affected by molecular absorption. Moreover, 
because of the pulsational mechanism, the effective temperature 
presents a strong variation in the luminosity cycle (as well as 
the associated black-body radiation). However, the emission fea- 
tures that we are analysing are often very intense (hence much 
higher than the local continuum). In the blue part of the spec- 
trum, where most of the Balmer lines are, the molecular absorp- 
tion is less, and the continuum change is assumed to be negli- 
gible with respect to the extent of the Balmer lines. However, 
we chose to use non-normalized data and a pseudo-continuum 
defined by ourselves. A polynomial fit of second order is cal- 
culated on the pseudo-continuum around the line and between 
-200km.s 'and 200km.s '. Then, we normalize the data with 
this fit. 

The polarization rates given in our study might reach unex- 
pectedly high values (several dozen percent). One could say that 
no polarizing mechanism may be able to polarize light as effi- 
ciently, but it has to be borne in mind that this rate is given as a 
percentage of the pseudo-continuum intensity, i.e. not as a per- 
centage of the total intensity (where the total intensity meant as 
the sum of the pseudo-continuum intensity and the intensity in 
the strong emission line), which is what is us ually given in theo- 
retic al works. This is the first time (except in iMcLean & Covna 
11978 ) that a polarimetric study is presented by analysing such 
intense emission. 

Compared to the work presented in lMcLean & Covn3l 19781 
and also devoted to omicron Ceti, our study is based on much 
more polarimetric observations collected at various phases. 
Since we present the evolution of the intensity and polarization 
with respect to the pseudo-continuum, those observations made 
through time are thus believed to be comparable. 

A threshold of detectability can be set at 10"^ (J. Morin, pri- 
vate communication). Also, a limiting constraint is the crosstalk 
between diff'erent Stokes parameters owing to instrumental prob- 
lems. It has been measured to be no more than 1 % (see NARVAL 
website). Knowing that, we chose 3% as a reasonable limit for 
polarimetric signatures to be considered as detections. 

2.3. Spectropolarimetric data around Balmer lines 

Figures |2] to |5] present the polarimetric data (in Hy and H6) 
of omicron Ceti for the best-sampled cycle of our survey (the 
third cycle). Our analysis is focussed on those two Balmer lines, 
which are the least affected by TiO absorptions. In Figs |2] and 
m the spectra of the four Stokes parameters I, Q, U, and V are 
displayed. The abscissa is taken as velocity in the ste llar rest 
fram e (omicron Ceti's radial velocity: V«=63.8km.s~ MWilsonl 
1 19631) . and the spectra were convolved with a Gaussian curve 
of FWHM = 4 km.s 'to slightly reduce noise and improve the 
quality of the data. In Figs [3] and |5] we also show u(q), which 
can be read as a polar graph with the linear polarization rate pii„ 
as the radial value and 20 as the angle. 

2.3.1. Hy 

From Figs |2] and [3] devoted to the data around the Balmer 
line Hy, the emission appears at i/j=2.95, hence slightly before 
the maximum light. During the full duration of our monitoring, 
the Hy emission line increases steadily and displays some well 
known features (Joy 1947, Gillet 1988, Castelaz et al. 2000). Its 
time evolution can be related to the presence of a shock wave 
emerging from the photosphere just before the maximum light 
and starting to propagate outward. While it first encounters a 



medium of decreasing density in the lower part of the atmo- 
sphere, the shock thus accelerates, before being progressively 
dimmed in the highest atmospheric layers because of the energy 
lost through radiation. While it is doing so, the column densi- 
ties of absorbants (located above the shock's front) decreases 
slightly. Those molecular absorptions are shaping the emission 
line to leave three peaks that can be easily identified from ^=2.95 
to (^=3.28. The most redshifted peak is constant (in intensity) 
through time, in contrast to the two others. The central peak 
evolves quite rapidly around the maximum light and then re- 
mains constant in intensity until the end of our observations. The 
bluest peak develops progressively and reaches (^=3.28 at an in- 
tensity that is five times higher than the one displayed at (p-2.95. 
On its blue edge, a shoulder is also present and disappears around 
(^=3.21. Each of these peaks is globally redshifted with time with 
hardly any change in the width. Throughout our observations, 
the total width of this emission is conserved through time at 90 
km.s"'. 

When the intensity starts to develop, the Stokes parameters 
q and u both present strong signatures. Stokes q has a positive 
value at phase (^=2.89, 3.05. It becomes zero and then negative 
through time, reaching its full amplitude at (f-3.l^. All the sig- 
natures in q spread over a velocity width equivalent to the one 
of the intensity emission and the features often present a line re- 
versal at their centre. The Stokes u signature remains positive 
throughout our monitoring. It reaches its maximum amplitude at 
^=3.09 and remains quite constant in shape and in amplitude un- 
til 1/3=3. 12. After the maximum light, the signatures on u present 
a strong change in shape, because strongly asymmetric and be- 
cause clearly becoming narrower with time. In a less definite 
way. Stokes v signatures are very weak and present structure 
variations through time. While the radiation is clearly getting 
linearly polarized during the shock's propagation (from i^=3.05 
to (/?=3.28), circular polarization is hardly detected except at the 
very end of our monitoring. 

The variation in the angle of linear polarization is quite obvi- 
ous when looking at the u(q) plots (Figure[3]l. At (f-3.Q5, the first 
phase of definite polarization, 6 is around 20°. Afterwards the di- 
rection of linear polarization rotates counterclockwise to reach 
90° at ^=3.28. After the maximum light, at phase ^=3.05, the 
linear polarization rate reaches 0.2 and stays around this value 
except for the last observation at i^=3.28, where it decreases to 
0. 1 . However, noticeable is that, when the intensity emission is 
still growing, the linear polarization rate decreases after its max- 
imum. Whatever is polarizing light works less efficiently at this 
specific moment. 

2.3.2. H6 

From Figs |4] and |5] which displays the data around the 
Balmer line H6, the intensity emission appears at (/?=2.95 similar 
to what is observed for Hy. It starts to develop rapidly after the 
maximum light and continuously increases during the monitor- 
ing to reach a maximum amplitude from ^=3.21 to ^=3.28. The 
maximum intensity of H6 is much higher than for Hy. However, 
the temporal evolution of H6 is different in the sense that the 
intensity stays constant during the last two phases in contrast to 
what happens in Hy. This is probably because H6 photons are 
the most energetic. Since the shock is weaker at that time, they 
should be the ones whose production should start to be reduced 
first, and H6 is also less mutilated by molecular absorption than 
the Hy line. One strong absorption at zero velocity severely af- 
fects the profile, while another one is barely noticeable on the 
blue wing around maximum light. The total width of the emis- 
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sion is constant through time at 75 km.s" , which is smaller than 
the one in Hy. For this line, the width and position are barely 
modified through time. 

As in Hy, from the maximum light up to the end of our ob- 
servations, signatures are clearly detected in Stokes u, q, but also 
in V. Stokes q varies from being slightly positive before the max- 
imum in luminosity, to nearly cancellation around (p-2.95, and 
then to negative values from (f-3.05. Then it starts to develop 
and is at a maximum from (^=3.18 until the end of our survey. 
The q signatures become narrow progressively and even present 
a small peak appearing at some phases on the bluest side. Stokes 
u is undergoing dramatic changes. After being roughly zero until 
(f-3.05, it develops positively, getting as high as 0.3 at (p=3.09. 
Then this positive peak recedes while a negative structure devel- 
ops, and finally at ^=3.28 u signature is fully negative. Over a 
large portion of our monitoring. Stokes u presents an "inverse 
P Cygni" structure, positive in the bluest part. This part of the 
spectrum corresponds to the one in intensity where the emission 
is the brightest. At last, the Stokes v parameter presents a very 
strong feature at phases (^=3.12 and 3.18, which almost disap- 
pears suddenly at following phases. 

The evolution of the linear polarization angle in the u(q) 
plane is slightly similar to Hy. For phases such as (/?=3.05 and 
1^=3. 12, the angle cannot be determined easily. However, we can 
definitely state that evolves from 45° at ^=3.09 to eventually 
settle down at approximately 100° at the last three observations 
(from (^=3.18 to ^=3.28). This is interesting when compared to 
Hy where does not settle down. The linear polarization rate 
is at first as low as 0.1 at (p-2.^9 and clearer at (p-3.Q5. It in- 
creases up to 0.5 at (/?=3.18 and 3.21 and then decreases to 0.4. 
This reduction in polarization seems to occur during the shock's 
weakening. 

During those phases. Stokes v shows up before settles 
down to some value. We know that the decrease in the inten- 
sity emission is d ue to the deceleration of the shock wave. In 
iGillet et alJ (fT98l . it is shown that the emission line in Ha in- 
creases until 1^=0.40 (the reference maximum is of course dif- 
ferent than ours) and then goes down. The change in the polar- 
ization behaviour described above thus certainly happens at the 
moment of this deceleration. 

2.3.3. Comparisons to previous cycles 

In Figures|6]to|9l we present some phases from the previous, 
but less well-sampled cycles. When we compare the same phases 
of different cycles, the intensities of Hy and H5 are roughly the 
same, in amplitude, width, and shape. Therefore, we can safely 
state that, for the considered phases of each cycle, the shock 
wave will produce intensity emissions of similar strength. 

In contrast, the u, q, and v Stokes parameters can strongly 
vary from one cycle to another Both symmetries and ampli- 
tudes of q, u, and v signatures are found to be altered in 
many cases, inducing a change in the linear polarization angle. 
Interestingly, the linear polarization rate seems to be roughly 
conserved through cycles. This allows us to stress the non- 
repeatability of the polarizing processes in action. Indeed, we 
know already that each cycle of the luminosity variation of omi- 
cron Ceti is not an exact copy of the previous ones (see AAVSO 
website). Considering that the shock wave is the consequence of 
the pulsation of the star, as the luminosity variation is, then it 
is no wonder that the polarizing mechanism linked to the shock 
wave does not behave exactly the same as before. However, as 
we noticed for the third cycle, the signatures' widths do not 
change from one cycle to another for a given Balmer line. 



3. Interpretation 

The appearance of an emission line indicates that the specific 
intensity of the back-lighting radiation (creating the continuum) 
is less than the source function inside the shock. This is linked to 
the intrinsic photon emission behind the shock's front. Emission 
lines can thus help in determining the properties of this area. 
Polarization also happens to be much stronger in the Balmer 
lines than in the pseudo-continuum around them. Its temporal 
evolution follows the evolution of intensity emission. It is there- 
fore believed that the polarizing mechanism is also intrinsic to 
the shock wave. 

In Figures |2] to |5] the evolution of intensity emission lines 
H6 and Hy is quite similar, leading us to think that the shock 
emerges from the photosphere just before the maximum light 
and is getting stronger toward the end of the third cycle of our 
survey. For both Balmer lines, the emission's width and central 
velocity (expressed in stellar rest frame) do not vary with time as 
previously said. This would imply that the shock is not acceler- 
ating much. Moreover, if the width of the emission remains the 
same, it probably means that the mechanism behind line broad- 
ening (micro-turbulence?) does not evolve with time. Willso^ 
(1976) already addressed the characterization of shock wave 
properties through emission lines width in Miras. In the case 
of an optically thin medium and when broadening line mech- 
anisms are not dominant, the width is an indicator of the flow 
velocity where the radiation is emitted. When the gas is opti- 
cally thick, we cannot give anything else but an upper limit for 
this velocity. Moreover, the mo del of shock wave we ch ose also 
influences this relationship. In iFadevev & GiUeti (l200l . shock 
waves models are used to create synthetic Balmer lines with dif- 
ferent values of temperature, density and upstream velocity. The 
models show a line doubling that is particularly obvious for Ha 
and just slightly obvious for H/3 while Hy stays single-peaked. 
On the other hand, the widths of the emissions vary clearly with 
density and upstream velocity. 

A few shock-related mechanisms of polarization are prob- 
ably operating. First of all, we consider a photospheric origin. 
The shock starts deep in the atmosphere, probably even be- 
low the photosphere. The photospheri c origin h as already been 
thought of for polarization (e.g. Harrington 1969) but never with 
a shock wave crossing it. Besides this, it is known that stars 
w ith low gravity, s uch as Mira stars, have large convective cells. 
In lChiavassa et al.l ( 1200 9), 3D models of red supergiants (RSG) 
photosphere were compared to interferometric observations in 
the near infrared. They conclude that RSG have large-scale gran- 
ulations undergoing a temporal vari ation. This kind of wor k 
has also been done for AGB stars dFrevtag & Hofned l2008l) . 
Thus, this convective structure could transfer its asymmetry to 
the propagating shock wave. Since the convective cells are very 
large, the asymmetry is not cancelled by integration on the stel- 
lar disk, and the shock wave may be aspherical. However, other 
sources of asymmetry ( e.g. high-layer asymmetries) might be of 
interest in this respect (Ragland et al. 2008, Pluzhnik et al. 200' 
Chiavassa et al. 2010, Wittkowski et al. 201 1). Fadeyev & Gille. 
( 1200 ih have shown that the compression rate r] behind a radiative 
shock can be written as 



:4-h3- 



F- 



+ 3- 



pUE, 



(1) 



where E, is the specific energy in the translational degrees of 
freedom, Ei„ the specific energy of excitation and ionization of 
atoms, and Fr the radiation flux. The subscript 1 refers to the 
values in the state ahead of the shock. Thus in the presence of 
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radiation, the gas compression ratio can reach values of tens or 
even hundreds, significantly greater than the compression ratio 
without radiation (rj = 4). With such high compression ratios, 
the level of turbule nce of the gas must be strongly amplified (see 
iGillet etal.iri998bl) . When the shock wave passes through the 
convection zone located around the photosphere, the convective 
cells are thus sharply compressed and heated, from the point of 
view of the shock, the convective structure is almost static. It is 
likely that the convective field is disrupted and even destroyed 
temporarily during the passage of the wave. Consequently, if the 
heterogeneous distribution of convective cells causes polariza- 
tion, then it should also be observed in the continuum and be 
disturbed during the photospheric passage of shock. It seems 
that the quality and temporal resolution of current observations 
are not yet sufficient to confirm this point directly. New observa- 
tions are needed. 

More local phenomena can be called upon. When the Mach 
number of shock wave increases, its structure quickly becomes 
complex. In particular, the charge separation within the shock 
front region may induce an electric field. The nature of the re- 
sultant interaction of the electric field with the gas flow has 
been investigated for a steady plan e shock in fully and par- 
tially ionized gases (see for instance lJaffrin|[T965[ iLu & Huang! 
[T974irVidar et al.|[T995h . The electric field is parallel to the di- 
rection of shock propagation. However, the presence of back- 
ground density or temperature gradients, for instance, causes an 
electrical current, inducing a curl of magnetic field. Thus the 
shockfront becomes a source of magnetic field, which pervades 
the unperturbed gas and the wake of the shock. These studies 
were limited to the region of viscous shock and did not consider 
the plasma state appearing later in the wake when the coupling 
with the radiation field can be fundamental. As discussed above, 
the solenoidal state of the postshock gas strongly increases in 
the presence o f a radiative field (shock turbulence amplification, 
see lGillet et a l. 1998b). Consequently, because the shock propa- 
gates in a non-uniform background plasma, a spontaneous gener- 
ation of a magnetic field is possible within the postshock plasma. 
From the induction equation, the equation describing the devel- 
opment of magnetic field is 



dB 
dt 



(2) 



The three terms on the right-hand side describe the convection, 
diffusion, and generation of a magnetic field. The self-generation 
magnetic field could be due to several local mechanisms such 
as ion-electron separation at the front or density and tempera- 
ture gradients in the plasma. Thus the dominant physical pro- 
cess must be determined to write the source term S. The mag- 
netic field is perpendicular to the shock propagation direction. 
Consequently, a linear polarization could be expected, although 
because of the complexity of the environment, a minor circular 
polarization could also be observed. 

The observation of the polarization of emission li nes ha s 
been a long- standing issue in sola r flare physics (Haud T98T1). 
in tok amak (iFuiimoto et al.|[T99 ^, gasdischarge ( Wald eri et alj 
Il999l) . and laser produced plasmas (Yonedaetal. 199?!). The 
emission polarization can come from different causes induced 
by the non-uniform background plasma. In particular, it has been 
known that anisotropic collisions of electrons with atoms pro- 
duce polarized radiations in plasma in which the collisional ex- 
citation and ionization processes dominate (impact polarization). 
Thus, the generation of polarized emission can occur during the 
recombination processes. In fact, there are few theoretical and 
experimental works about this matter, and devoted to radiation 



fronts, laboratory ionizing shock tubes, or radiative shock waves. 
In conclusion, we expect that the observed polarization in emis- 
sion lines of Mira could have been caused by the spatial non- 
uniformity in the electron distribution function that was induced 
by the self-generated magnetic field created in the shockfront. It 
is not possible within the framework of this study to definitively 
explore the physical origin of this polarization. Only a quantita- 
tive approach may allow us to answer 

Tables |2] and [3] sum up the evolution of spectropolarimetric 
values for Hy and H6 for the 3rd cycle. Considering only I pro- 
files, it is difficult to state at which moment of the cycle the shock 
reaches its maximum intensity because the emissions do not un- 
dergo the same molecular absorption. Then, polarization could 
help determine this moment. If we consider the local mecha- 
nism, the turbulence would be at its maximum during phases of 
most intense emissions, inducing the strongest polarization. It is 
interesting to see that linear polarization appears during a longer 
time than the circular one. Stokes V is rarer and mostly present 
when the linear polarization is quite strong, i.e. at about phase 
0.20. Indeed, this is consistent because when the shock becomes 
strong, the turbulence amplification is at a maximum, and con- 
sequently the hydrodynamic motions are more and more tridi- 
mensional. It is interesting to note that phases around 0.20 are 
quite common in atmospheric models for merging shock fronts 
Clreland et al..,2011i) . 

Table 2: Time evolution of spectropolarimetric values in Hy. 
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Table 3: Time evolution of spectropolarimetric values in H6 (see 
Table|2]i. 
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To know which of those two phenomena are likely to occur, 
it is necessary to make similar observations of other types of 
radially pulsating stars that have no convection, but are known 
to host an atmospheric shock wave. The detection of polariza- 
tion in the Balmer emission lines of such stars would then give 
credit to the local hypothesis. On the other hand, if the polar- 
ization arises from the convective asymmetry, it would also act 
on the continuum of the photosphere's radiation. At that point, 
the global hypothesis would prevail. A temporal survey of the 
polarization of this continuum would then be necessary to see 
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whether it varies Hke the Balmer lines polarization presented in 
TablesElandEl 



4. Conclusion 

We performed the first spectropolarimetric monitoring of the 
Mira star omicron Ceti, focussing on Balmer emission lines. 
Thanks to these observations, we were able to detect variable 
linear and circular polarization signatures associated to those hy- 
drogen lines. The temporal variation of those signatures is cor- 
related with the evolution of the intensity of the emission lines 
known to be formed in the radiative wake of a shock wave prop- 
agating throughout omicron Ceti's atmosphere. 

Several potential polarizing mechanisms are considered. 
AGE stars host very large photospheric convective cells induc- 
ing overhall asymmetry. Since the shock is formed near the 
photosphere, it is possible for it to interact with the convec- 
tive structure and thus to gain the photospheric asymmetry. 
Interferometric observations at minimum light, i.e. when the 
shock is high in the atmosphere, would help to compare the 
photosphere's and the shock's geometrical features. On the other 
hand, our local hypothesis states that a charge separation induced 
by the shock's passage along with the turbulence believed to ex- 
ist behind the front of the shock could produce a local magnetic 
field through which the radiation from the wake becomes polar- 
ized. 

Unlike omicron Ceti, other Mira stars can be rich in carbon 
or have no companion. They would thus be interesting to observe 
to determine how general this evolution of the line polarization 
is. This study might also be extended to other kinds of pulsating 
variables stars, such as RV Tauri or RR Lyrae, known to host 
radiative shock waves propagating throughout the atmosphere. 
If pulsating stars with no convection (such as y6 CMa stars) were 
still showing such line polarization, it would favour the local 
hypothesis. 
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Fig. 2: Hydrogen Balmer lines signatures of the third cycle of our monitoring of omicron Ceti in Hy (full Stokes). Successive 
signatures are artificially offset vertically for legibility. Phases are indicated on the right (see Table[TJ. 
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Fig. 3: Hydrogen Balmer lines signatures of the third cycle of our monitoring of omicron Ceti in Hy (u(q)). Phases are indicated in 
the topleft of each frame (see Table[T]i. 
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Fig. 4: Hydrogen Balmer lines signatures of the third cycle of our monitoring of omicron Ceti in H6 (full Stokes). See Figure|2l 
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Fig. 5; Hydrogen Balmer lines signatures of the third cycle of our monitoring of omicron Ceti in H6 (u(q)). See Figure[3] 
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Fig. 6: Hydrogen Balmer lines signatures of cycles 1 (thick lines) and 2 (thin lines) of our monitoring of omicron Ceti in Hy (full 
Stokes). See Figure|2] 
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Fig. 7: Hydrogen Balmer lines signatures of cycles 1 and 2 of our monitoring of omicron Ceti in Hy (u(q)). See Figure[3] 
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Fig. 8: Hydrogen Balmer lines signatures of cycles 1 (thick lines) and 2 (thin lines) of our monitoring of omicron Ceti in H5 (full 
Stokes). See Figure|2] 
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Fig. 9: Hydrogen Balmer lines signatures of cycles 1 and 2 of our monitoring of omicron Ceti in H6 (u(q)). See Figure[3j 



